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1. Introduction 

Methods based upon coherent states [1] have proved to be widely useful in investigating 
the ground state properties of algebraic models. In algebraic models, the Hamiltonian 
and all physical operators are constructed from the elements of a Lie algebra, usually 
from a bosonic realization of XJ(n). Such models have been applied extensively to 
the spectroscopy of many-body systems, including nuclei [2] and molecules [3]. The 
coherent states of an algebraic model are obtained by repeated action of a general 
linear combination of boson creation operators on the vacuum state. As variational 
trial states, the coherent states allow the estimation of the ground state energies and 
properties, yielding results which become exact in the infinite boson number limit [4, 5]. 
They are also essential in defining the classical limit for the model [5], providing the 
geometric coordinates or dynamical variables of the model. 

Coherent states may also be used to study the intrinsic excitation modes of an 
algebraic system [6-15], through the construction of coherent states orthogonal to the 
ground state. However, application of this method has generally been limited to excited 
states involving only one intrinsic excitation quantum, or at most two [15], due to 
the complexity of calculating the necessary matrix elements. In the present work, a 
general expression is obtained for the matrix element of an m-body operator between 
coherent states constructed from multiple orthogonal linear combinations of boson 
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creation operators. This allows the coherent state formalism to be applied to states 
possessing an arbitrarily large number of intrinsic excitation quanta. 

The necessary definitions for the basic ground state (condensate) coherent state 
are first presented, and the results for expectation values with respect to this state are 
reviewed (section 2). The general result for excited coherent states is then established 
(section 3). As a simple illustration, calculations are carried out for excited states in the 
SO(3) dynamical symmetry limit of the molecular two-dimensional vibron model [16] 
(section 4). 



2. Condensate coherent state 



An algebraic model based upon a bosonic realization of U(n) is obtained by defining 
n bosonic creation operators b\, b 2 ,. ■ ., b^ n obeying the canonical commutation relations 
[bi,bj] = 0, [bi,bj] = Sij, and [b\,bj] = 0. The set of all possible bilinears of a creation 
and an annihilation operator then forms a basis b\bi, b\b 2 , . . ., b^bn for the algebra. The 
physical operators of the model, such as the Hamiltonian, the angular momentum, and 
transition operators, are constructed as polynomials in the b\bj. These operators act on 
the states of the Fock space created by b\, b 2 , ■ ■ •, b\ t . 

The condensate coherent state is defined in terms of the "condensate boson" 
creation operator, which is a general linear combination 

B\ = a x b\ + a 2 b\ + ■■■ + a n b\, (1) 

with complex coefficients satisfying the normalization convention Yll=i a *i a i = 1- The 
normalized condensate coherent state is then 

\N;a 1 ,...,a n ) = ^=(Bt) N \0). (2) 



This state is an eigenstate of the total number operator A^ = Y17=i K^i- 

The expectation value of a one-body or two-body operator with respect to the 
condensate (2) was deduced by Van Isacker and Chen [17], using arguments based upon 
formal differentiation. Here let us derive an explicit result for the expectation value of 
an arbitrary m-body operator, since the results of section 3 can be obtained as a natural 
extension. First, for an annihilation operator b r , note that [b r , =a r , from which the 
relation 

[b r ,(Bt) N } = Na r (Bt) N - 1 (3) 
follows by the product rule for commutators. The action of b r on the condensate is thus 

b r \N;a!, . . . ,a n ) = VNa r \N - 1; cci, . . . , a n ). (4) 
Repeated application yields the expectation value of an arbitrary m-body operator, 

mm m 

(N; a n \ (J[ b\?j (J[ b r ^j\N; a 1 ,...,a n ) = N^Y[a* r ,a ri , (5) 

i=l i=l i=l 

where the underlined superscript indicates the falling factorial [18], defined by m- = 
m(m — 1) • ■ • (m — n + 1). 
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3. General coherent state 

For the study of excited states, it is necessary to consider coherent states which are 
orthogonal to the condensate state. These are constructed using multiple different 
coherent boson species B\ (s — 1, . . . , S), defined as linear combinations 

B\ = a S)1 b{ + a Sj2 b\ H 1- ®s,nbl (6) 

of the basic boson creation operators. (The procedure for obtaining values of the 
coefficients a s ^ appropriate to a given model is discussed in, e.g., ref. [11].) The 
coefficients a S;i are chosen to obey the orthonormalization convention Yli=i a t>,i a s,i — 
S S ' :S . Consequently, the coherent bosons satisfy canonical commutation relations 
[B s ,, B s ] = 0, [B s ,, Bl] = Sj,,, and [B\,, B\\ = 0. The different B\ represent the ground 
state condensate boson and one or more orthogonal excitation modes. The normalized 
multi-species coherent state is 

\N 1 ...N s ) = (fl-^(Bl) N ')\Q). (7) 



s=l v * 

The coherent state is an eigenstate of the total number operator, of eigenvalue N = 
J2s=i-N s - States of different coherent boson occupation numbers N ± , . . ., N s are 
orthogonal. 

The matrix element of an arbitrary m-body operator with respect to two multi- 
species coherent states can now be deduced following the approach of section 2. The 
commutation relations generalize to [b r , B^ s ] = a s ^ r and 
s s s 

[br, 11(^1 = E N «*r [U^sf^] ■ ( 8 ) 

s=l t=l s=l 

The action of b r on the multi-species coherent state is thus 

s 

KIN, ...N S ) = J2 y/Ntat^Ni - S ltt ) ■■■(N s - 6 s , t )). (9) 
t=\ 

This is analogous to the result of equation (4), but with a separate term arising from 
the action of b r on each species of coherent boson contributing to the coherent state. It 
is useful to define a counting function v(ti, . . . , t m ; s) = YlT=i ^U,s, giving the number of 
the ti which are equal to s. Then, for a product of annihilation operators acting on the 
coherent state, 

m S S I to 

(n*.)w •••"»>= e [nv^^Kn ««.«) 



i=l ti,...,t m =l s=l i=l 

x | (iVi - v(h, . . . , t m ; 1)) • • • (Ns - v(h, ...,t m ; S))}. 



(10) 
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Table 1. The number of terms in general contributing to the sum in equation (11), 
not vanishing due to the delta symbol constraint, for various values of S and m. In 
practice, some of these terms may also vanish as a consequence of zero values for the 
falling factorials or a St i coefficients. 



s 




Contributing terms 




1-body 


2-body 


3-body 


4-body 


2 


2 


6 


20 


70 


3 


3 


15 


93 


639 


4 


4 


28 


256 


2716 


5 


5 


45 


545 


7885 



The matrix element of an arbitrary m-body operator (m > 1) with respect to two 
arbitrary multi-species coherent states is the inner product of two such expressions, 

m m 

(K---N' s \(jJb^(J[b rt )\N 1 ...N s ) 

i=l i=l 

S S i ; m (11) 

= e [n 5 N's-^-»* v n> s-^] (n a u^) , 

t , -^,...,t f 7n _ y s 1 i 1 

where the abbreviations u' s = v(t\, . . . , t' m ; s) and v s = u{t\, . . . , t m ; s) have been used. 

Three stages are involved in evaluating the matrix element of a general operator: 
reexpression of the operator in terms of normal-ordered m-body terms, evaluation of 
the matrix elements of these by equation (11), and simplification of the result. For 
complicated operators or if many coherent boson species are involved, these steps can 
most effectively be carried out though computer-based symbolic manipulation. A few 
useful special cases of equation (11) are summarized in the appendix. 

The multiple sum in equation (11) nominally contains S 2m terms. However, for 
fixed numerical values of the N s and N' s , many of the terms vanish identically due to the 
restriction on indices imposed by the product of Kronecker delta symbols. A summary 
of the number of nonvanishing terms for various S and m is given in table 1. A given 
species s of coherent boson is overannihilated when N s — v{s) < or N' s — u'(s) < 0, 
yielding a vanishing falling factorial in equation (11). Thus, additional terms vanish if 
the expression is evaluated for a small value (< m) of any of the N s or N' s , as typically 
occurs when the lowest-lying excited states are considered. If the multiple sum in 
equation (11) is instead to be evaluated with the N s and N' s retained as variables, all 
terms involving the same product of falling factorials may be collected. This product is 
identical for terms with the same values of all the z/(s) and v'(s). Since < u(s) < m 
and YL 8 =i u ( s ) = m i an d similarly for z/(s), the number of distinct terms after collection 
is the square of the number of possible partitions of m over S bins. 

Multi-species coherent states of the form (7) are also encountered as the condensate 
states of systems involving multiple constituents, each separately conserved. An example 



Application of the coherent state formalism to multiply excited states 



5 



from nuclear physics is the proton-neutron interacting boson model (IBM-2) [2], in which 
proton pairs (created by s\ , d\ 2 i d\ ^ d\ , d\ +l , d\ +2 ) and neutron pairs (created by 
sl , d\_ 2 i dl_ ± , dl , dl +l , dl +2 ) are separately conserved. The physical operators are 
constructed from the elements of a Lie algebra Ui(n)®U2(n)®- • •, and a condensate state 
with good boson number for each constituent is constructed as oc (b1 1 ) Ni (b1 2 ) N2 ■ ■ ■ |0). 
Since the condensate bosons B\ p (p = 1, 2, . . .) are constructed from disjoint sets of boson 
operators, the expectation value of an m-body operator in general factorizes into the 
product of simple expectation values of the type (5) (e.g., ref. [19, (CI)]). However, the 
general result (11) for the multi-species coherent state matrix element can provide the 
simplest framework for computer-based symbolic evaluation [20]. 



4. Intrinsic excitations of the two-dimensional vibron model 

The two-dimensional vibron model [16] is the U(3) algebraic model, describing a system 
containing a dipole degree of freedom constrained to planar motion. The basic example 
of such a system is a triatomic linear bender molecule, but the model is easily extended 
to more complex molecular systems. The U(3) algebra is realized in terms of the three 
bosonic operators o~\ rj, and rj, which satisfy canonical commutation relations. It is 
convenient to define circular bosons rj_ = ^(rjiirj) / \[2 [21-23], such that the operators 
o~\ T+, and t_ carry 0, +1, and —1 units of two-dimensional angular momentum. The 
physical operators include the angular momentum operator Z = t|t + — t1t_, the dipole 
operators D± = ±v^2(t"1o" — crV T ), and the quadrupole operators Q± = \[2t±t^. 
The U(3) algebra contains the subalgebra chains [16] 

, x f U(2) D SO(2) , x 

U 3 D < V ; 12 

\S0(3) D SO(2). 

The dynamical symmetry associated with the U(2) chain yields spectra matching those 
of the cylindrical oscillator (Poschl- Teller potential), while the dynamical symmetry 
associated with the SO (3) chain yields spectra like those of the displaced cylindrical 
oscillator (Morse potential). The SO (3) limit is used in the following illustrations, as 
the less trivial case. The SO(3) subalgebra is spanned by {£>+, D-, 1} and has quadratic 
Casimir operator W 2 = (D+D- + D-D + )/2 + P. The SO(2) subalgebra is simply the 
two-dimensional angular momentum algebra, containing /. The simplest Hamiltonian 
with SO(3) dynamical symmetry is H = —W, which has eigenvalues 

E(N, v, I) = -N(N + 1) + 4v[(N + 1/2) - v], (13) 

with v = 0, 1, ... , LiV/2j and I = —(N — 2v), —(N — 2v) + l,..., +(N - 2v) [figure 1(a)]. 

The condensate boson and an orthogonal excitation boson for the two-dimensional 
vibron model may be defined as 

Bl(r) = ^=s(J + rrl) B\{r) = -^=={-ra^ + rj). (14) 
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V V 

Figure 1. Comparison of excited coherent state estimates of observables (curves) with 
exact values (circles) for the SO (3) dynamical symmetry of the two-dimensional vibron 
model, with N = 100. (a) Schematic SO (3) level energy diagram, with arrows indicating 
the dipole (solid) and quadrupole (dashed) transitions considered in panels (c,d). 
(b) Energies of all SO (3) representations, (c) Intensity of the (I = 1) — ► (I = 0) 
dipole transition within each representation, (d) Intensity of the (I = 2) — ► (I = 0) 
quadrupole transition within each representation (solid curve, filled circles) and of the 
(I = 2) — > (I = 0) and (Z = 0) — (Z = 2) Av = — 1 transitions between representations 
(dashed curve, open circles). An l 2 energy splitting is included in panel (a), to provide 
visual separation of levels within an SO (3) representation. All observables are plotted 
rescaled by 1/N 2 . 



(See ref. [11] for further discussion of the choice of boson operators in the vibron model.) 
The general excited coherent state is 

\NN x ;r) = 1 =[Bt(r)] N - N *[Bt(r)] N *\0). (15) 

' ' y/(N - N X )\N X V K J 

The expectation value of the SO (3) Casimir operator with respect to an arbitrary excited 
coherent state is evaluated using equation (A. 2), yielding 

(NN x ;r\W 2 \NN x ;r) 

4 (16) 
= 2[N + N X (N - N x )\ + - ^ 2)2 [N(N - 1) - QN X (N - N x )]r 2 . 

The equilibrium value of r is found by minimization of the variational energy 
(NO; r\H\N0; r), giving r — 1. With this value of r, the excited coherent state energies 
are 

E(N,N X ) = -N(N + 1)+4N X [N-N X ]. (17) 



Application of the coherent state formalism to multiply excited states 



7 



If the intrinsic excitation number N x is identified with the SO(3) quantum number v, 
comparison of equations (13) and (17) shows that these expressions differ only by a term 
of order 1/JV. Thus, to leading order in 1/N, the coherent state estimate reproduces 
the excitation energies for all excited states, as illustrated in figure 1(b). 

Transition strengths may be estimated using coherent states, from the squared 
matrix element of the transition operator, I (NN' x \T\NN x ) 2 . The coherent state 
\NN X ) is not an angular momentum eigenstate, so the resulting estimate for the 
transition strength between two intrinsic excitations is effectively averaged over the 
many states of different angular momenta constituting that excitation. This coherent 
state estimate thus cannot be expected to provide the exact transition intensity between 
any two particular angular momentum eigenstates. It does, however, indicate the general 
magnitude of the transition strengths and the overall dependence on excitation quantum 
number, and it can be quantitatively accurate if the angular momentum dependence of 
transition strengths is weak. (Alternatively, angular momentum eigenstates may be 
projected from the coherent states [9,11,24], but this requires additional machinery 
beyond simple evaluation of an m-body operator matrix element.) 

Dipole infrared transitions in the two-dimensional vibron model are, to leading 
order, induced by the operators D±, and quadrupole Raman transitions are induced by 
the operators Q± [3, 16]. The strengths of transitions within an intrinsic excitation are 
estimated from the expectation values 

(NN X ; r\ D ± \NN X ; r) = -2(N - 2N X )^-^ 

/MM in \MM \ 1 N * + ( N ~ 
(NN X ; r\ Q± \NN X ; r) = --^ ^— , 

and those between succesive intrinsic excitations from 



(18) 



y»2 



1 + r 2 



(N(N X - 1); r\ D± \NN X ; r) = -^/( N -N X + 1)N X 
(N(N X - l);rj Q± \NN X ; r) = ~ y/{N - N x + 1) jV^ 



(19) 



v/2 v v x ' ' x l+r 2 ' 

The coherent state estimates of transition intensities for the SO (3) dynamical symmetry 
are obtained from these equations with r = 1. The estimates are plotted for iV = 100 
in figure l(c,d), together with exact values obtained by numerical diagonalization, as 
functions of the excitation quantum number v (or N x ). The coherent state estimate for 
dipole transitions closely reproduces the strength of the angluar momentum 1 — > 
transition within an SO (3) representation [figure 1(c)]. Dipole transitions between 
different representations are forbidden, and the coherent state estimate indeed vanishes. 
Quadrupole transition strengths exhibit greater angular momentum dependence within 
a representation, and the coherent state estimate is consequently less accurate. The 
strengths of quadrupole transitions involving the low angular momentum members of 
the representations are reproduced to within ~ 5%, except at the highest intrinsic 
excitation quantum numbers [figure 1(d)]. Note that the angular momentum 2-^0 



Application of the coherent state formalism to multiply excited states 



8 



and 0-^2 transitions between two representations (Av = ±1) differ in strength, and the 
coherent state estimate consistently behaves as their average. 

5. Conclusion 

The present results serve as a basis for application of the coherent state formalism 
to states with an arbitrary number of intrinsic excitation quanta. This process yields 
estimates of eigenvalues and operator matrix elements for excited states valid to leading 
order in 1/N. The illustration provided was to a dynamical symmetry limit of a 
simple model, but the coherent state analysis will likely be most useful when applied to 
transitional Hamiltonians, for which analytic results are not otherwise available. 

The coherent state formalism has in the past provided not only a quantitative 
calculational tool but, perhaps more importantly, a method for obtaining qualitative 
understanding of the equilibrium properties and fundamental modes of a system. Most, 
if not all, of the raw numerical results of the coherent state formalism can also be 
obtained by numerical diagonalization. It is thus the latter, interpretational aspects of 
the coherent state formalism, and the explicit analytic forms obtained for the parameter 
dependences of observables, which have proved most useful. The present results for 
multiply excited states thus might most productively be used in investigating the general 
nature of the evolution of a system's properties with excitation energy. 
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Appendix. Special cases of the general matrix element 

Some useful special cases of the general matrix element (11) are given here. If the 
coherent states involve only two species of coherent boson (S = 2), then the 1-body and 
two-body operator matrix elements are 

((Aq - 1)(N 2 + 1)| 0% |AqA 2 ) = y/N^ + 1) a^a^ 

(AqA" 2 | |AqA 2 ) = Aq al r ,ai tr + N 2 a 2r ,a 2tr (A.l) 

((Aq + 1)(N 2 - 1)| blb r |AqA 2 ) = v/(Aq + l)^ 2 oqVa 2ir 



(A.2) 
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{{Nl _ 2 )(N 2 + 2)| bl,bl,b ri b r2 |iViiV 2 > = V / W^«; ir ,a;,v«i 1 n«i,r 2 
((N, - 1)(N 2 + 1)| b\,b\,b rx b r2 \N X N 2 ) 

= (Ni - l)y/N 1 (N 2 + 1) (al r2 ,al ri , + a\ r2 ,a*^ ril )a 1>ri a 1>r2 
+N 2 ^/N l (N 2 + 1) a* 2r ^a* 2ri ,{a 2ri a lr2 + a lin o; 2)J . 2 ) 
(A^iA^ 2 | b\,b\,b rx b r2 \N X N 2 ) = a\ r2 ,a\ ri ,a lri OL lr2 + JV 2 - Ci* 2r2 iCi* 2ri ,a 2ri a 2r2 

+N X N 2 (al r . 2 ,a 2 , ri , + a^al^a^a^ + a 2 ^a^ r2 ) 
{{N\ + 1)(N 2 - l)\bl,bl,b ri b r2 \N!N 2 ) 

= iViv^Vi + l)N 2 a* lr2 ,a\ ril (a 2ri a lr2 + a lri a 2r2 ) 

+ (N 2 - 1)^(^1 + l)-^2 (a*l,r 2 ' a 2, ri > + ^.ra'^l.nO^.n^.ra 
((TVi + 2)(iV 2 - 2)| bl,bl,b ri b r2 \N,N 2 ) = y/iN, + WN^<,r 2 '< ri ^2,r^2,r 2 , 

with all others zero. For an expectation value (all N' s = N S ), equation (11) simplifies to 

m m 

(N 1 ---N s \( 1 [[bl^(j[b ri )\N 1 ...N s ) 

i=i i=i 

s s to (A. 3) 

References 

[1] Zhang W, Feng D H and Gilmore R 1990 Rev. Mod. Phys. 62 867 

[2] Iachello F and Arima A 1987 The Interacting Boson Model (Cambridge: Cambridge University 
Press) 

[3] Iachello F and Levine R D 1995 Algebraic Theory of Molecules (Oxford: Oxford University Press) 

[4] Gilmore R, Bowdcn C M and Narducci L M 1975 Phys. Rev. A 12 1019 

[5] Gilmore R 1979 J. Math. Phys. 20 891 

[6] Bohr A and Mottelson B R 1982 Physica Scripta 25 28 

[7] Bijkcr R and Diepcrink A E L 1982 Phys. Rev. C 26 2688 

[8] Leviatan A 1985 Z. Phys. A 321 467 

[9] Leviatan A 1987 Ann. Phys. (N. Y.) 179 201 
[10] Kuyucak S and Morrison I 1988 Ann. Phys. (N. Y.) 181 79 
[11] Leviatan A and Kirson M W 1988 Ann. Phys. (N.Y.) 188 142 
[12] Kuyucak S and Morrison I 1989 Ann. Phys. (N. Y.) 195 126 
[13] Leviatan A and Kirson M W 1990 Ann. Phys. (N. Y.) 201 13 
[14] Leviatan A and Ginocchio J N 1991 Phys. Lett. B 267 7 

[15] Garci'a-Ramos J E, Alonso C E, Arias J M, Van Isackcr P and Vitturi A 1998 Nucl. Phys. A 637 
529 

[16] Iachello F and Oss S 1996 J. Chem. Phys. 104 6956 

[17] Van Isacker P and Chen Jin-Quan 1981 Phys. Rev. C 24 684 

[18] Graham R L, Knuth D E and Patashnik O 1994 Concrete Mathematics: A Foundation for 

Computer Science 2nd ed (Reading, MA: Addison- Wesley) 
[19] Ginocchio J N and Leviatan A 1992 Ann. Phys. (N. Y.) 216 152 



Application of the coherent state formalism to multiply excited states 



10 



[20] Caprio M A and Iachello F in press Ann. Phys. (N. Y.) 
[21] Perez-Bernal F private communication 

[22] Perez-Bernal F, Santos L F, Vaccaro P H and Iachello F submitted Chem. Phys. Lett. 

[23] The definitions r| = T(Vj ± ir+)/v^ [21,22] are used, rather than r| = (rj ± ir^/v^ [16]. 
This choice is necessary for _D± to be the physical dipole operators and provides a closer 
correspondence with the treatment of the 3-dimensional vibron model [11]. 

[24] Pcicrls R E and Yoccoz J 1957 Proc. Phys. Soc. (London) A 70 381 



